ITP Budapest Rep. 537 
BI-TP 97/45 
November, 1997 



Coupled Gap Equations for the 
Screening Masses in Hot SU(N) Gauge Theory 

A. Patkos 1 * P. Petreczky 1 '^, Zs. Szep 1 * 
x Dept. of Atomic Physics, Eotvos University, 

H-1088 Puskin 5-7, Budapest Hungary 
2 Fakultat fiir Physik, Universitat Bielefeld, 

P.O. Box 100131, D-33501 Bielefeld, Germany 



Abstract 

Coupled 1-loop gap equations are studied numerically for non-Abelian electric and 
magnetic screening in various versions of the three-dimensional effective gauge models. 
Corrections due to higher dimensional and non-local operators are assessed quantita- 
tively. Comparison with numerical Monte-Carlo investigations suggests that quanti- 
tative understanding beyond the qualitative features can be achieved only by going 
beyond the present treatment. 
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1 Introduction 



An outstanding consequence of heating non-Abelian gauge fields is the screening of static 
chromo-electric and -magnetic fields. Electric (Debye) screening is generated by non-zero 
Matsubara modes both for Abelian and non-Abelian gauge fields and its value at leading 
order has been known since a long time. The dynamics of the non-Abelian zero modes is quite 
complex [|l| and it leads, at least in the Schwinger- Dyson approach to the static magnetic 
gluon two-point function, to the generation of a magnetic screening mass Qj. Non-zero 
magnetic mass was obtained from all variants of the magnetic gap equation in the effective 
three-dimensional gauged Higgs and also in pure gauge theories 0, f| §], H f| . 

The spectra of screening masses in the SU(2) Higgs model were also studied in lat- 
tice Monte-Carlo simulations by measuring various correlation functions of gauge invariant 
operators [J71 || as well as the gauge boson propagator in fixed (Landau) gauge JTU|] . Elec- 



tric and magnetic screening are simultanously accessible to four-dimensional Monte-Carlo 



simulations [O. The magnetic mass obtained from the gauge boson propagator in Landau- 
gauge is rather close to that obtained from the gap equation, while the gauge invariant 
gluonic correlation functions yield several times larger mass. A consistent interpretation of 
the situation emerging in the weakly coupled Higgs+Gauge model was suggested in (see 
also |l3j) in the framework of the so-called constituent model, where the mass measured 
in fixed gauge correlators corresponds to the constituent magnetic gluon mass of the con- 
fined three-dimensional theory. The idea of the constituent gluon offers a natural effective 
degree of freedom, which dominates the static part of the free energy of the gluon plasma 



14| . It should be noticed in this context that the relation between gap equation and the 
resummation of the free energy has been recently discussed in [JOj, [RJ |HJ for the scalar field 
theory. 

For pure hot SU(N) theory the situation is more complicated since for the physically 
interesting range of the temperature the value of the gauge coupling g is close to unity and 
it is not clear whether the standard argument for dimensional reduction applies. 

The aim of the present paper is to investigate how far the screening masses can be 
described by means of the technique of gap equations. Since the hiearchy of scales 2nT >> 
gT » g 2 T fails to hold in pure SU(N) theory for the temperature range of interest we 
should investigate a coupled set of gap equations for all the screened modes and determine 
the corresponding screening lengths simultanously. 

The most straightforward way to do this would be to derive gap equations in the full four- 
dimensional theory. However, it is not known how to generalize the by-now well-established 
three-dimensional gauge invariant resummation || [(J to four dimensions. There are gauge 
transformations which might mix static and non-static modes, therefore the resummation of 
the static modes only, which was suggested in \TE\ violates gauge invariance. Since screening 
masses are static quantities it is natural to calculate them in the framework of an effective 
three-dimensional theory which is, however, valid only up to scales k ~ gT. It was shown 
that the accuracy of the description of such theories is improved if in the action of the 
effective theory beside superrenormalizable operators one also takes into account higher 
dimensional and non-local operators [[19] . When deriving the gap equation in the framework 
of three-dimensional adjoint Higgs model it is important to address the question whether the 
symmetric or the broken phase is the physical one. In Refs. |7|, |2(J it has been argued that 
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symmetric phase is the physical one, the study of the Debye screening also suggests that the 
presence of the A condensate is physically unfavourable, therefore in what follows we will 
assume no A backround. 

In this paper we shall assume only the separation of the static and non-static scales in 
the pure SU(N) gauge theory, that is we assume that g(T) « 2tt. In this way we will 
be interested in the derivation and the solution of the coupled gap equations of the electric 
and magnetic static fields. To our best knowledge our paper represents the first attempt for 
this simple generalisation of the gap-equation approach, where till now electric and magnetic 
screening masses have been discussed separately. The closest to the spirit of our investigation 
is the analysis of the sensitivity of higher order corrections of the electric screening mass to 
the existence and the value of a magnetic screening scale by Rebhan |H|] . 

The paper proceeds as follows: in section 2 we shall calculate screening masses from the 
local superrenormalizable effective theory (i.e the adjoint Higgs model) of the hot SU(N) 
gauge theory using a gauge- dependent resummation scheme (namely, i?^-gauge). The tem- 
perature and gauge-parameter dependence of the results are discussed in detail. Gauge 
invariant resummation schemes are used for the derivation of similar coupled gap-equations 
in section 3. In section 4 we shall analyze the effect of higher dimensional and non-local 
operators in the gap equation. In section 5 detailed comparison with the results of earlier 
and the most recent numerical investigations will be presented together with our conclusions. 



2 Gauge non-invariant resummation scheme 

This scheme for the evaluation of the magnetic mass was first suggested in || and for the 
Debye mass in ||21|| . It was noticed in Ref. M that magnetic mass obtained in this scheme is 
gauge dependent and therefore cannot be regarded as physically meaningfull. However, even 
in gauge invariant resummation schemes the value of the magnetic mass cannot be defined 



unambigously, because it depends on the specific resummation scheme [22]. Therefore it 
seems interesting to calculate the magnetic mass in a gauge non-invariant scheme just to 
compare the amount of gauge dependence with the ambiguity of gauge invariant approaches. 

The three-dimensional effective Lagrangian relevant for 1-loop calculations can be written 
with i?£ gauge fixing as 

L = \F*F* + \ (( D t A ) 2 + mlA-Al) + 
\m\A\Al + (5 4 c a Ac a + m 2 G c a c a ) + j-^Ai) 2 + L ct , 
Let = \{m 2 m - ml)A a A a - \m\A\A* - m 2 G c a c a (1) 

with 

DiAl = diA% - gsf^A^Al, 
F§ = diA] - djAl + gzf abc A\A ). 

Here we have added and substracted a mass term for Ai, Aq and the ghost fields with their 
exact screening masses. m|) is the tree-level (from the point of view of the effective theory) 
Debye mass for A , which was generated during the procedure of the dimensional reduction. 
The other parameter of the effective theory is the three-dimensional gauge coupling ^3 which 
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appears in the definitions of the covariant derivatives and the field strength tensor. At 1-loop 
level one has m 2 m = g 2 NT 2 /3 and g\ = g 2 T, where g = g(T) is the gauge coupling of the 
original theory. The propagators can be read from the quadratic part of the Lagrangian 
and can be found in the Appendix, where one also finds some details of the evaluation of 
the relevant Feynman diagrams contributing to the different 2-point functions. It should be 
noticed when performing the resummation of the pure gauge sector with the unique mass 
term of mass tjit the longitudinal and transverse gluons acquire different masses, which are, 
however related by m^ = y/^rriT {ttit is the transverse and is the longitudinal mass). 
It is also possible to perform the resummation by introducing independent masses for the 
longitudinal and transverse gluons, but then the corresponding gap equations will have only 
complex solutions. The gauge boson self-energy can be decomposed as 

Hij{k) = (6ij - --^-)n T (A;, m T , m D , m G ) + -^-U L (k, m T , m Dl m G ). (2) 

In the Appendix we give the expression of in terms of a few fundamental three-dimensio- 
nal loop-integrals. These integrals are easily evaluated and an explicit but very cumbersome 
functional form can be written for the longitudinal and transversal projections of the polar- 
isation matrix. 

The self energy for A was first calculated in [plfl : 



HooO, m D , m T ) = m 2 m + 



2 , g 2 h 



Itt 



. 2{m 2 -k 2 -m 2 /2) 

-mD — mx H — - — a — r U — L arctan ■ 



(3) 



k mo+m-T 

+ (k 2 + m 2 ) ( (arctan ^ arctan * ) + - 1) J- 

V U '\ m r k V m D + ^/lrn T m D +m T ' VVS I mT 

The on-shell gap equations now can be written as 

m 2 ^ = Uxik = imT,mT,mD,mG), 

m 2 D = n 00 (A; = im D ,m D ,m T ), (4) 
m 2 G = U G (k = im G ,m T ,m G ). 

On the mass-shell n o(& = imp, mo, tut) is gauge parameter independent, but Urik = 
im T ,m T ,m D ,m G ) and U G (k = im G ,m T ,m G ) do depend on the gauge fixing parameter, 
therefore the masses obtained from this coupled set of gap equations are gauge dependent. 

In the following numerical investigations we shall consider the case of the SU(2) gauge 
group. The 4-dimensional coupling constant is taken at scale ft = 2ttT, where ft is the MS 
scale and 1-loop relation for the gauge parameter of the effective theory is used. To set the 
temperature scale we use the relation T c /A^r = 1.06 obtained from numerical simulation of 
the finite temperature SU{2) gauge theory [|lT| . 

The temperature dependence of mp in a specific gauge is plotted in Fig.l. As one can see 
from the plot mjy receives 30% positive correction compared to the leading order result, while 
the magnetic mass stays very close to the value calculated by Buchmuller and Philipsen ||, 
given below in Eq.@. Since the masses are gauge dependent in this approach, it is important 
to investigate the dependence of the screening masses on the gauge parameter £. It turns 
out that one gets real values for the masses from the gap equations only if £ € [1, 5) . The 
dependence of the screening masses on £ in this range is shown in Fig. 2. One can see, that 
the ^-dependence of my in this interval is 40 %, while for mo it remains in the 10 % range. 
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Figure 1: The temperature dependence of m D /m D0 (solid line) and m T /m pp for Feynman 
(£ = 1) gauge, moo is the leading order result for the Debye mass and m pp is the value of 
the magnetic mass obtained by Buchmuller and Philipsen for pure SU (2) gauge theory. 
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Figure 2: The dependence of m D /m D0 (solid line) and m T /m p on gauge parameter £ at 
T = 10 4 T C , m D0 is the leading order result for the Debye mass and m pp is the value of the 
magnetic mass obtained by Buchmuller and Philipsen for pure SU (2) gauge theory. 
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3 Gauge Invariant Approach 



Gauge invariant approaches for the magnetic mass generation in three-dimensional pure 
SU (N) gauge theory were suggested by Buchmuller and Philipsen (BP) || and by Alexanian 
and Nair (AN) ||. The approach of AN uses the hard thermal loop inspired effective action 
for the resummation of the magnetic sector. The approach of BP using a gauged u-model, 
goes over to the SU(N) gauge theory in the limit of infinitely heavy scalar field. Till now 
only these two gauge invariant schemes are known to provide real values for the magnetic 
|22 . 



mass 



In these approaches the gauge boson self-energy is automatically transverse and there 
is no need to project the transverse part from the polarisation tensor. The corresponding 
expression for the on-shell self-energy reads 



H T (k = im T , m T ) = Crriq 



(5) 



where 



C 



(6) 



'^[f ln3-lUAr 
^[|ln3-i],SP. 

Since we are interested in calculating the screening masses in the three-dimensional 
SU(N) adjoint Higgs model, U T (k,m T ) should be supplemented by the corresponding con- 
tribution coming from A fields. This contribution is calculated from diagrams d) and e) of 
the Appendix to be 



47T 



+ 



k 2 + Am 2 D 
4k 



arctan 




(7) 



It is transverse and gauge parameter independent, it also does not depend on the specific 
resummation scheme applied to the magnetostatic sector. It should be also noticed that it 
starts to contribute to the gap equation at 0(g 5 ) level in the weak coupling regime, thus 
preserving the magnetic mass scale to be of order g 2 T. This is the reason why no "hierarchy" 
problem arises in this case, at least for moderate g values. 

The self energy of A depends on the specific resummation scheme. For BP resummation 
it reads 



Il 00 (A;,m D ,m T ) 



m, 



DO 



+ 



9 2 N 
4tt 



. 2(m 2 D -fc 2 - m 2/2) 

—tjid — tht H — - 12 — z — LJ — L arctan 



rriT+rriD 



(k 2 +m 2 D ) 



— ( — m T + \ arctan - 



This expression is different from the expression of n o calculated in the gauge non-invariant 
approach (see eq. ) but its analytic properties and on-shell value is the same as of (|3]). 
For the resummation scheme of AN the self-energy expression of Aq coincides with ([|) if it 
is evaluated at £ = 1. The coupled set of gap equations now can be written as 



m\ = CrriT + SH A °(k = imT,rri£)), 
m 2 D = U 00 {k = im D ,m D ,m T ). 



(9) 
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Figure 3: The temperature dependence of the scaled Debye mass for BP resummation scheme 
(solid) and for the AN resummation scheme (dashed). The scaling factor is m D0 . 




Figure 4: The temperature dependence of the scaled magnetic mass for BP resummation 
scheme (solid) and for the AN resummation scheme (dashed). The scaling factors are m^ p 
and mj> N , respectively. 
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The temperature dependence of m D obtained from this coupled set of gap equations is shown 
in Fig. 3 for both schemes, where we have again normalized the Debye mass by the leading 
order result, m^. The temperature dependence of the magnetic mass is shown in Fig. 4, 
where we have normalized ttit by the value of the magnetic mass obtained for pure three- 
dimensional SU(2) theory, in the BP (AN) gauge invariant calculations |5], [J. As one can 
see the contribution of A Q to the magnetic mass is between 1 and 3%. From Figures 3 and 
4 it is also seen that the temperature dependence of the screening masses is very similar to 
the temperature dependence of the respective leading order results. 



4 Contribution of non-Local Operators to the Gap 
Equation 

In the previous section the gap equations were derived for an effective local superrenor- 
malizable theory. In this case the effect of non-static modes in the 2-point function were 
represented by the thermal mass for Ao and by the field renormalization factors which relate 
3d fields to the corresponding 4d ones. It was shown in [^, |24| that when one performs the 
procedure of the dimensional reduction in gauges the parameters of the superrenormaliz- 
able effective theory are gauge independent, only the expressions of 3d fields in terms of 4d 
ones depend on the gauge parameter. However, this does not hold for higher dimensional 
and non-local operators, which are generally gauge dependent. At 1-loop level the only 
diagrams contributing non-locally to the gap equation are those which have two non-static 
line inside the loop, diagrams with one static and one non-static line inside the loop are 
forbiden because of 4-momentum conservation. Therefore at 1-loop level the non-locality 
scale (27tT) _1 is much smaller than the relevant length scales. 

In general the non-static contribution to the static 2-point function n Atl/ (/c = 0, k) can 
be written as 

An-(/c = o, k) = 5°5°Mk) + fiistfa - ^)n(fc), 

^oo(k) = m 2 m + ai (fi, i)k 2 + T 2 E^ =2 a„(0(2&) 2n , 



IT»(fc) = b 1 ^,Ok 2 + T*En=2UO(£rY^ ( 10 ) 

where the coefficents a n and h n can be calculated for arbitrary n. The first two terms in the 
expressions of n o and first term of U(k) are already included into the 3d effective theory as 
part of the tree level mass and the definition of the 3d fields in term of 4d fields. The last two 
sums will contribute to the 3d effective lagrangian as quadratic non-local operators. There 
are also higher dimensional operators as well as non-local 3- and 4-point verticies in the 
effective lagrangian, however, since we restrict our interest to 1-loop gap equations these are 
not important for us. Their contribution would correspond to 2 or higher loop contribution 
in the full four-dimensional theory. 

We have estimated by direct numerical evaluation of the infinite sums in ([TO) the contri- 
bution of non-local operators to be less than 1% in the temperature range T = (3 — 10 4 ) T c , 
thus neither their contribution nor their gauge dependence is essential. 



2n 
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5 Conclusion 



In the present paper we have made an attempt to extract the electric and magnetic 
screening masses from the coupled set of gap equations of the three-dimensional SU(N) ad- 
joint Higgs model considered as an effective theory of QCD. The screening masses have been 
studied using gauge non-invariant as well as gauge invariant resummation schemes. In the 
gauge non-invariant formalism we have observed rather strong gauge parameter dependence, 
therefore the results extracted from it are not very informative. It is still interesting to note 
that in Feynman gauge (£ = 1) the results for the magnetic mass are rather close to those 
obtained from the gauge invariant resummation scheme of Buchmuller and Philipsen. In 
gauge invariant treatments we have compared two different resummation schemes, that of 
Buchmuller and Philipsen (BP) and one proposed by Alexanian and Nair (AN). Qualitativ- 
ely these two resummations lead to similar results, but in the BP scheme one has smaller 
magnetic mass and larger Debye mass than in the AN scheme. 

Let us summarize our view on the interaction of the electric A and the magnetic Aj 
fields. In both schemes one can see that the dynamics of Aq is largely influenced by the 
magnetic sector, however, no similar feedback on the magnetic sector is seen, the magnetic 
masses calculated from the coupled gap equations provide screening masses which are 1% 
smaller than evaluated in the pure gauge theory. This fact suggests that the role of the 
adjoint scalar field is similar to that of the fundamental Higgs field, because the magnetic 
mass calculated in the symmetric phase of SU(2) Higgs theory using lattice Monte-Carlo 
simulation with Landau gauge fixing is also roughly the same as in pure gauge theory (jTIJ . 

Finally we compare our results with recent Monte-Carlo data for the screening masses 
obtained in 4d finite temperature SU(2) gauge theory |TT]. The data on the magnetic mass 
found from this simulation in the temperature range T = (10 — 10 4 )T C can be fitted well by 
the formula % = 0.456(6)g 2 T. This value of the magnetic mass is considerably larger than 
what one obtains from the magnetic gap equation and 3d simulations, where the results are 
approximately rriT = 0.28g 2 T for the BP scheme, 0.38g 2 T for AN scheme and 0.35g 2 T for 
3d simulation. Also our results are rather close to these values, as we have demonstrated in 
Sections 2 and 3. 

The data from Monte-Carlo simulation for the Debye mass in the above mentioned tem- 



perature interval can be fitted using the following leading order-like anzatz J 1.69 (2) g(T)T 



lH , which means that the Debye mass is roughly 1.6m do, where m^o is the leading order 
result. The gap equations at the same time, as one can see from Fig. 3 give (1.2 — 1.3)mno, 
depending on the resummation scheme. While there is no quantitative agreement between 
masses measured in Monte-Carlo simulation and those obtained from the gap equation, the 
temperature dependence of these masses in the temperature interval T = (10 — 10 4 )T C seems 
to follow the temperature dependence of the leading order result. 

At the moment we have no explanation for the reasons of the discrepancy between the 
results of the 4d simulation and the 3d gap equations and simulation. The observations may 
imply that either higher loop contributions to the gap equation are important and therefore 
the dynamics of A is not faithfully reflected by 1-loop gap equations or the dimensional 
reduction is not valid. In |25| the self-consistency of the gap equation has been examined 
in the entire momentum space in gauge invariant way. This investigation reveals some 
inconsistencies of the 1-loop gap equation, which are conjectured to be removed by higher 
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order calculation. It remains an open question how far these inconsistencies influence the 
pole mass. We believe, however, that further studies in this direction, which include 2-loop 
gap equations |26| and lattice study of the effective adjoint Higgs model [^7] will clarify part 
of these problems. 
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Appendix 

The propagators of different fields can be read off the quadratic part of the Lagrangian. 
These exact propagators are listed below. The gauge boson propagator: 



kikj \ 1 k{kj £ 

k 2 j k 2 + k 2 k 2 + m 2 L ' 



A#) = U, " H ) j^—r + ^T^-r- (H) 



where mt = y/^mT, the propagator for the adjoint scalar field A : 

1 



D Ao (k) = TT—^i (12) 



m 



2 • 

D 



and finally, the ghost propagator: 



m = tt—^- (13) 



k 2 + 



rn, 



G 



The Feynman diagrams contributing to the 2-point functions of the relevant fields are 
shown below. In these diagrams every line corresponds to a resummed propagator. The wavy 
line corresponds to the gauge bosons, the solid one to A and the dashed one to the ghost. 
The corresponding analytic expressions can be written in terms of the following standard 
integrals: 



I(k,r, s) 



d d p 1 
(2?r) d (p + k) 2r p 2s 

k d-2(r+ S ) r(r + s -|) r(f- s )r( 



I m (k,r,s) 



( 47r )d/2 r(r)r(s) r(d-s-r) 

d d p Pm 

(2n) d (p + k) 2r p 



(14) 



,2s 



k d-2(r+ S ) r (r + s - f) r (I + 1 - s) r ( 



(4 n )d/2 T(r)T(s) r(d + l-s-r) 



(15) 
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(27r) d (p + k) 2r p 2s 

k d-2(r+s) [k 2 r (r + s - 1 - 1) r (f + 1 - s) r (f + 1 - r) 



J(k,m 1 ,m 2 ) 
J m (k,m 1 ,m 2 ) 



( 47r )d/2 



<i 3 p 



r(r)r( a ) 



T(d + 2 -s-r) 



r r+ s -f r f +2- 



s .r(|-r 



r(r)r( a ) 
1 



T(d + 2- s-r) 



(2tt) 3 ((p+ A;) 2 + m 2 )(p 2 + m 2 ) 8tt' 
(2tt) 3 ((p + /c) 2 + m 2 ) (p 2 + m 2 .) 



J mn (k,m 1 ,m 2 ) = J 



8tt 
d 3 p 



mi — m 2 k 2 



m\ — wij 



k 2 



2k 2 



A 



(2%) 3 ((p + k) 2 + ml) (p 2 + m 2 2 ) 



,2\2 



8tt 

k m k n 
8tt 



mi ^m 2 , + (A; 2 + m 2 — m|) 
~2~ + 8P A 



2^2 



mi \k 2 m\ + 3 (A; 2 + m 2 — m 2 ) 
2P 



8A; 4 



1 (mi — m 2 ) {k 2 + m 2 — m 2 , 



8tt 



4A; 2 



/c 2 



where A — \ arctan ( — j — ) 

k \mi+m 2 / 

j(m) 
/ mn (m) 



d 3 p 



(2ir) 3 p 2 + m 2 

d 3 p PmPn 



f an 

J (2n 



m 

'4tt' 



m 2 j(m). 



(27r)V + m 2 3 

These integrals were evaluated using dimensional regularization. 

The diagrams contributing to the 2-point functions of Ai are the following: 





a 
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with the analytical contribution: 



n«7(*) 



mn 
a6 



J mn (A;, m G , m G ) + k m J n (k, m G , m G ) 



n (c) 7(A:) = -\g*N5 ab {5, 



4 2 

- (m ^ + 2 fc2) ( J(k, m T , 0) + J(fc, 0, m T )) 

Tflrp 



(23) 
(24) 



2k 2 {Am 2 T + k 2 ) jn . / fc 2 \ .. . 

+ 5 -J(k,m T ,m T ) -2 ( 1 + — j(m T ) 



~\~k m k n 
+ 
+ 



(m 2 T + k 2 ) 2 . . ( k A 
4 J(fc,m r ,0) + I 



A' 2 



A: 1 



4m 4 
3 A; 2 



4mf, 



6 ) J(k, ttit, mr) 
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4 + ^ ^ + ^ ) Ak, 0, m T ) + ^jK) + X^X^, 1, 1) 
£■ A- £- 4 / £- 4 P 



4mn 



4m^ 



m 2 j(mr) + 4^ WA; ' 1 ' 1)+ l4^ + 4 m 2 



8 J mn (k,m T ,m T ) 



rn 



2 ( k 2 \ 2 

—lmn(m T ) -IH 2" (Jmn(k, 0, m T ) + J m „(fc, ra T , 0)) 



+ 



m 5 , 

—^-3{m T ) + ^kralnik, 1, 1) - ( -^y + 1 ) ( + 3 ) k m J n (k,m T , 0) 



( k A \ ( k A 

+ I 1 I A; m J n (fc, 0, m T ) + I — + 4 



my 
A; 2 



rrin 



rrin 



rrin 



k m Jn 



(k,m T ,m T ) | 



-g 2 NS c 



4 



A- A- 



(ml + k 2 ) 



m, 



(J(k, 0, m T ) - J(fc, mi, m T )) 



(ml + m? + A; 2 ) 



4m^m 2 



+ 



^SA^m^ + 7m 4 -. + 2m^m 2 — (m 2 + A; 2 ) ^ J(/c, mi, m T ) 
+m 2 L j{m T ) — 7m^j(m L ) + (A; 4 — 6Aj 2 m|, — 7m^ J(k, 0, m r ) 
-A; 4 /(A;, 1, 1) + (m\ + A; 2 ) 2 J(k, m L , 0) 

m 2 T k m k n j{m L ) + A; 4 ( J mn (fc, m L , 0) - I mn (k, 1, 1)) 

+ (A; 2 + m|) 2 (J mn (A;,0,m T ) - J mn (k,m L ,m T )) + m 2 T l mn (m L ) 

L ^-k m J n (k,m L ,m T ) ^-^k m I n (k, 1, 1) 



m\m\ 



2 2 



m T m 



2 \ 2 ~"~ I ^rn^n\i^) m Li^>) i 2 2~~ 



A- A- 



m 2 L 



j{m L ) + 



mj,m" L 



k m J n (k,0,m T ) 
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-2m 2 L] (m L ) + k 4 (I(k, 1, 1) + J(k, m L , m L )) 
(k 2 + mf) 2 J{k,m L , 0) - (fc 2 - m|) 2 J(Jfe, 0, m L ) 



ml 



I mn {k, 1, 1) - J m „(fc, m L , 0) - J mn (fc, 0, m L ) + J m „(fc, m L , m L ) 



k 2 ( k 2 \, /7 , k 2 ( k 2 \, 

T ~T + 1 k m J n (k, m L , 0) 2" ~ 1 k rnJn{k, 0, 771l) 

mi Vmi y m^ v^i ; 

fc 4 

(M,i))| 

mi 



n«7(*) 

n (e) 



u6 



4 J mn (fc, m D) m D ) + 4k m J n (k, m D , m D ) + k m k n J(k, m D , m D ) 



mn 
ab 



g 2 NS ab S mn j(m D ). 



(25) 

, (26) 
(27) 



Diagrams a), b) and c) were calculated in |28j using the Landau gauge (£ = 0). The results 
of these calculations coincide with ours if in the above formulas one sets £ = 0. The diagrams 
contributing to the 2-point function of A are given by diagrams f) and g): 



_ _ 




g 



The corresponding analytical contribution is: 



n^I(fc) 



-g 2 N8 ab 8V 
1 



]{m T ) 



k 2 + m 2 T 



m 



3 



(m ^ + 2 fc2)2 J(A:,m g ,0) (28) 
mm 



r 



2k 2 (m? + m 2 D ^j + k 4 + [m 2 ^ — m 2 ^ j J(fc, m£>, my) — j(rrio) 



g 2 N5 ab 5^5 0v \ \(m 2 L + m 2 D + k 2 ) j(m L ) + (m 2 , + A; 2 ) 2 J(fc, m D , 0) 



1 

m' L 

m D + ^(^j m D, m L ) 



nto) 



/if 
«6 



g jyoab 



2j(m T ) +U(m L ) 



(29) 



Finally the single diagram contributing to the ghost 2-point function is given by diagram h): 

S ~ \ 



h 



12 



The corresponding analytical contribution is: 

1 f r™2 (h.2 i ,^2 

Z ab (k) = --g 2 N5 ab ^ 



ml ~ + mh) j(m T ) - ^ t ^ J(k, m G , 0) 



+ 



'(k 2 



2 \ 2 



4mt 



m 



G) 



k 2 — m% m% 



1 



4 J(k,m G ,m T ) - ^j(m G ) 



{k 2 + trig) 2 k 2 + m; 



Am 2 L 



2 m 2 \ 1 
^ + J ( fc ' m G' m ^) - 4^(^g) 



2 



4m| 



-J(k,m G ,0) 



3m 2 , 



k 2 + m 



G 



4m| 



(30) 
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